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ABSTRACT 

We show that the relativistic wind of the Crab pulsar, which is commonly thought 
to be invisible in the region upstream of the termination shock at r < rs ~ 0.1 pc, 
in fact could be directly observed through its inverse Compton (IC) 7-ray emission. 
This radiation is caused by illumination of the wind by low-frequency photons emitted 
by the pulsar, and consists of two, pulsed and unpulsed, components associated with 
the nonthermal (pulsed) and thermal (unpulscd) low energy radiation of the pulsar, 
respectively. These two components of 7-radiation have distinct spectral characteris- 
tics, which depend essentially on the site of formation of the kinetic-energy-dominated 
wind, as well as on the Lorentz factor and the geometry of propagation of the wind. 
Thus, the search for such specific radiation components in the spectrum of the Crab 
Nebula can provide unique information about the unshocked pulsar wind that is not 
accessible at other wavelengths. In particular, we show that the comparison of the 
calculated flux of the unpulsed IC emission with the measured 7-ray flux of the Crab 
Nebula excludes the possibility of formation of a kinetic-energy-dominated wind within 
5 light cylinder radii of the pulsar, R w > 5i?L- The analysis of the pulsed IC emission, 
calculated under reasonable assumptions concerning the production site and angular 
distribution of the optical pulsed radiation, yields even tighter restrictions, namely 
Rw ^ 30i?L- 
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1 INTRODUCTION 

The Crab pulsar is a powerful nonthermal machine, acceler- 
ating plasma in the form of a relativistic wind that carries 
off most of the rotational energy of the pulsar. 

At a distance of about r = rs ~ 0.1 pc the wind is ter- 
minated by a standing reverse shock, which accelerates the 



electron s up Lu energies 10 eV, and landumizcs Lhci i pilch 



angles gees fe Qunn 107-1 r Kennel fc Coroniti 19 84 ). This 



results in formation of a bright synchrotron source in the 
region downstream of the shock. The synchrotron radiation 
of the Crab Nebula is well studied in a very broad frequency 
range, from radio to hard X-rays. Its general spectral and 
spatial characteristics are satisfactorily explained by the rel- 
ativistic ma gneto hydrodynamics (MHD) model of Kennel 
& Coroniti (1981). Remarkably, the latter provides also a 



reasonable explanation, even in its simplified (spherically 
symmetric) form, for the detected very high energy (VHE) 
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7-rays as a result of inverse Compton (IC) scattering of rela- 
tivistic electrons in the ambient low-frequency photon fields. 
This implies that the study of the TeV IC radiation of the 
Crab Nebula, combined with synchrotron X-ray emission, 
can yield unambiguous information about the relativistic 
electrons and the neb ular magnetic field in the downstream 
region of the shock J d^_Ja^c^&HaKhirg_^992 ^_ Stcpanian 



1995 ; Atoyan & Aharonian, 1996 



Hillas et al. 199^ )7 

Although very important, this information unfortu- 
nately does not tell much about the origin and character- 
istics of the wind, i.e. about the region between the pulsar 
magnetosphere and the shock. It is generally believed that 
this region, where almost the whole rotational energy of the 
pulsar is somehow released in the form of kinetic energy of 
the wind, cannot be directly observed. This has a simple 
explanation. Although the wind electrons may have an en- 
ergy as large as ~ 10 13 eV, they move together with the 
magnetic field and thus do not emit synchrotron radiation. 
This explains the f act that the region upstre am of the shock 
is underlurninous (Kennel & Coroniti 1984). However, this 



statement is valid only for the synchrotron radiation of the 
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wind. In fact, the wind could be directly observed through 
its IC radiation. Indeed, the IC 7-radiation of the wind elec- 
trons is unavoidable because of the illumination of the wind 
by external low-energy photons of different origin. There are 
three isotropic photon field populations that contribute ef- 
fectively to the production of IC 7-radiation of the nebula in 
the downstream region: the nonthermal (synchrotron) and 
thermal (dust) radiation of the Crab Nebula itself, and the 
2.7 K microwave background radiation (Atoyan & Aharo- 
nian 1996). Formally, all these photon fields could also serve 
as a target for the IC scattering of the wind before the shock. 
However, the fluxes of 7-rays emitted through these chan- 
nels appear to be well below the sensitivity threshold of 
7-ray telescopes. Meanwhile, in the pulsar vicinity, namely 
within approximately 100 light cylinders, the efficiency of 



produc jion ot UJ 7-rays dramatically increases because ot 
the existence of intense low energy radiation of the pulsar 



itself. The radiation consists ot two, pulsed and unpulsed 
components, the latter being modulated at the period of 
the pulsar. 

In this paper we show that at certain circumstances con- 
cerning the position of formation of the particle dominated 
wind, the geometry of the flow and the Lorentz- factor of the 
bulk motion of the wind, the fluxes of the IC 7-radiation 
of the wind could be sufficiently high to enable detection 
by present and forthcoming space-borne and ground-based 
7-ray telescopes. Moreover, the distinct spectral features of 
this radiation could allow effective separation of the "wind" 
component of radiation from the heavy background, which 
consists of unpulsed radiation of the Crab Nebula at very 
high (TeV) energies and pulsed radiation at low (GeV) en- 
ergies. We argue, that even upper limits obtained in such a 
study could provide unique information about the origin of 
the pulsar wind. 



2 CHARACTERISTICS OF THE WIND 

2.1 The total particle ejection rate of the wind 

The wind from the Crab pulsar carries away most of the 
energy of rotation of the pulsar. The energy released in the 
form of electromagnetic emission of the pulsar, which peaks 
at gamma-ray energie s, does not e xceed 1 per cent of the 
total rotational losses (Axons, 1996), therefore it can be ne- 
glected in the energy balance of the wind. 

The fluxes of the energy and the angular momentum of 
the wind consist of two parts. One of them corresponds to 
the matter and another corresponds to the electromagnetic 
field. The total flux of the kinetic energy of particles can be 
presented as 



Skin = Nmc < 7™ >, 



(1) 



while the flux of the angular momentum of the matter is 
equal to 



ikin = Nm < r^ w v v > . 



(2) 



Here N is the total rate of particle ejection from the pulsar 
magnetosphere, < 7 W > is the average Lorentz-factor of the 
wind, r is the distance to the axis of rotation v v is the com- 
ponent of the velocity of plasma propagation in the direction 
of rotation. Hereafter we assume that the wind consists of 



only electrons and positrons. The total rate of particle ejec- 
tion can be estimated within the models pair production 
in the pulsar magnetosphere: 



N„ 



(3) 



where S cap = 2nRlQ/c is the total area of the polar caps 
of the neutron star where the roots of the open magnetic 
field lines are placed, n± is the density of the particles in 
the primary beam, Q is the angular velocity of rotation of 
the pulsar, R, is the radius of the pulsar and the factor 
A takes into account the multiplication of particles because 
the development of electromagnetic cascades in the mag- 
netosphere. In the inner gap models, the electromagnetic 
cascade is initiated by a beam of e lectrons accelerated up to 



the Lorentz-facto r 7 gap ~ 2 ■ 10 7 (Ruderman & Sutherland 



1975; Arons 1985). The density of particles in the beam is of 



the order of G oldreich- Julian density n± = tigj (Goldreich 



Julian 196!:) determined as 



(OB) 

2-7rec 



(4) 



Owing to the electromagnetic cascade in the pulsar magne- 
tosphere, the number of particles increases by a factor of 
A and, correspondingly the Lorentz-factor of particles de- 
creases by the same factor. 

The calculations by D augherty & Harding (1982) and 



Gurevich & Istomin (1985) show that for the Crab pulsar 
this factor is of order of 10 4 . Note that these early calcu- 
lations took into account only multiplication resulting from 
the cascades supported by two processes - the curvature 
radiation of electrons and e pair production of 7-rays of 
these photons in the magnetic field. Meanwhile, the process 
of the Compton scattering of electrons on the soft thermal 
emission of the neutron star plays, most pro bably, a n on- 



negligible role in the cascade development ( Arons 1998 ), 
thus this effect should be taken into account in the esti- 
mates of A. In any case, the uncertainties in the model pa- 
rameters and assumptions does not allow an accurate the- 
oretical estimate of A, but instead give a broad range of 
possible values of A between 10 3 and 10 5 . Remarkably, a 
rather accurate estimate of < 7 W > in the upstream flow 
can be derived from the analysis of the nonthermal high 
energy radiation of the downstream region. Indeed, the ex- 
planation of the spectrum of synchrotron X-ray emission by 
the wind electrons accelerated (redistributed) by the ter- 
mination shock requires a power-law injection spectrum of 
the electrons Q(E) oc E~ a with a ~ 2.4 and a cutoff of 
the spectrum below E* ~ 150 — 200 GeV. Also, the inter- 
pretation of the X-ray and TeV 7-ray emissions within the 
synchro- Compton model of the Crab Nebula allows one to 
derive, with very good accuracy, the average magnetic field 
in the downstream region and the total luminosity in shock- 
accelerated electrons, B ~ 2-10~ 4 G, and W 



(for review see Aharonian & Atoyan (1998)). The obvious 



3xl0 38 erg/s 



conservation laws concerning both the number and the to- 
tal energy of the of relativistic electrons in the downstream 
and upstream regions (i.e. before and after the shock accel- 
eration) gives < 7 W >= §5± (E,/m c c 2 ) ~ 1.3 ■ 10 6 . The 
accuracy of this estimate depends on the possible range 
of variation of the parameters a and E* that still fit the 
data, and is estimated smaller than factor of 4. Corre- 
spondingly we find the injection rate of the wind particles 
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N = W I < 7w > m c 



2.8- 10 3 



This implies that for 



JVgap = 2 ■ 10 34 the A-factor should be close to 10 4 . 

Thus, the cascade multiplication of the primary elec- 
trons accelerated in the pulsar magnetosphere leads to the 
formation of an plasma with parameters 



N = 



and 



7w0 = 



\B D. 2 Ri 



7gap 

A ' 



The Lorentz factor of this initial wind is close to 7^,0 
while the kinetic energy flux of the wind is 



(•») 



(6) 



i(r 



E « -y wo mc 2 N (7) 

For the Crab pulsar this flux is estimated as 1.6 x 10 35 erg/s. 
It is much lower than the total rotational losse s of t he pul- 
sar E rot = 4 x 10 38 erg/s (see e.g. Shklovsky ( Jl968| )). This 
leads to a conclusion that the most part of pulsar's rota- 
tional ener gy is c arried off by the electromagnetic field (see 
e.g. Arons (1996)). The state of the wind is characterized by 
so-called <r-parameter determined as the ratio of the electro- 
magnetic energy flux to the kinetic energy flux of particles in 
the wind. At a w > 1 the wind is Poynting-flux dominated. 
At tj w < 1 the wind is kinetic-energy-dominated. Initially 
the wind ejected from the pulsar magnetosphere is Poynting 
flux dominated, because o" w o = E TO t/Eo ~ 2.5- 10 3 . The esti- 
mates of the wind parameters in the outer gap model gives 
essenti ally the same magn etization parameter (Cheung & 
Cheng |L994 |Coroniti 199q ). 

On the other hand, the explanation of the characteris- 
tics of the nonthermal radiation of the Crab Nebula requires 
that a w be between 10 -3 and 10~ 2 in the wind region up- 



stream of the ter mination shock (Rees & Gunn 1974 



Kennel 



& Coro liti 1984 ) . Thus the magnetization parameter a w de- 
creases by several orders of magnitude on the way from the 
light cylinder to the pre-shock region. Therefore it is difficult 
to avoid a conclusion that the wind is additionally acceler- 
ated in a some region beyond the light c ylinde r. 

The theory by Kennel & Coroniti ( 1984 ) is based on 
the assumption of ideal MHD flow of the plasma after the 
terminating shock wave. All dissipative processes are ne- 
glected, with exception of the cooling o f the plasma due 
to the sync hrotron emission. Lyubarskii (1992) and Begel- 
man (1998) argued that the theory of Kennel & Coroniti 
(1984) is likely not complete. The magnetic field after the 
terminating shock is mainly toroidal. It is very well known 
in p lasma physics t hat such configuration is strongly unsta- 
ble ( [Bateman 198C| ). The instabilities can essentially change 
the physics of the flow in the Nebula. The basic process is 
the fast dissipation of the toroidal magnetic field with con- 
version of its energy into the energy of the relativistic parti- 
cles. The instability drives the plasma towards equipartition 
of energy between the magnetic field and the matter. This 
process will inevitably be accompanied by acceleration of 
particles. The dissipative processes provide the dynamics of 
the plasma, in good agreement with the observed velocity 
of the expansion of the outer edge of the nebula, even for 
the parameter cr w ~ 1. We emphasize that, although very 
reasonable, this argument nevertheless does not solve the 
problem of the wind acceleration, because even in this case 
we have to transform wind with <r w o ~ 10 3 into wind with 



cr w ~ 1. Moreover, the wind with <j w ~ 1 has almost the 
same characteristics as we presented above. The only ex- 
ception is that the Lorentz factor of the wind with <j w ~ 1 



is twice as small as that of the wind with o\, 



10- 



In 



the limits of uncertainty of the multiplication factor A this 
difference is not important, however 

2.2 The energy spectrum of the wind electrons 

Here we assume that the wind with a <C 1 is formed in some 
'acceleration region' at a distance to the axis of rotation R w . 
We also assume that the plasma axially isotropically fills all 
the open field lines. Thus, the wind is not modulated in the 
azimuthal direction. 

For calculation of the spectra of emission we need in- 
formation about spatial and energetic distribution of par- 
ticles of the wind. Let us first summarize the information 
about these characteristics following directly from observa- 
tions. The average Lorentz-factor of the wind in the regime 
of a < 1 is determined as 



< 7w >= 



E r , 



N 



(8) 



with a typical value ~ 10 . 

The kinetic-energy-dominated wind is believed to be 
cold, because the region of the flow of this wind is ob served 
as an 'underluminous' region (see Kennel & Coroniti ( 1984 ) 
and referenced literature) . Otherwise, a hot wind would pro- 
duce remarkable synchrotron emission, which would contra- 
dict the existence of underluminous region. 

There is definite latitudinal inhomogeneity of the wind. 
The observations by ROSAT of the torus of X-ray emission 
in the Crab Nebula clearly de monstrate that ene rgy flux in 
the wind varies with latitude (Hester et al. 1995). This fact 
means that the density of the wind, its Lorentz-factor and 
the toroidal magnetic field should depend on the latitude. 

It is believed that the existence of an X-ray torus implies 
that the wind flows predominantly in narrow disc-lik e sector 
along the e quator with opening angle less than 30° (Hester 
ct al. 1995). However, one should take into account the fact 
that optical emission is uniformly distributed in the Nebula. 
This means that there should be a component of the wind 
at high latitudes. As the luminosities in optical and in X- 
rays are comparable the energetics of this high latitude wind 
component is comparable with the energetics of the equato- 
rial component which produces the X-ray torus. Moreover, 
the evidence that the wind at high latitudes e xists close to 



the p ulsar follows also from HST observations ( Hester et al. 
199S| ). In our calculations we assume that the particle flux 



in the wind is isotropic. However to be consistent with ob- 
servations in X-rays we assume that the Lorentz-factor and 
the toroidal magnetic field of the wind depend on latitude. 

There are no model independent estimates of the lat- 
itudinal dependence of the Lorentz-factor of the wind and 
toroidal magnetic field. We use here the characteristics of 
t he wind obtaine d in the model of an axisymmetric rotator 
(Bogovalov 1997) . These charact eristics are also valid for 
oblique rotators (Bogovalov 1999). 

Below we assume that 

7w = 7w0 + 7max COS 2 (q) , (9) 

where 7 max = o"wo7o is the maximum Lorentz-factor of the 



© 0000 RAS, MNRAS 000, 000-000 



4 S.V. Bogovalov & F.A. Aharonian 



plasma on the equator(a = 0); 70 ~ 10 3 (see equation ^ 
and a w o ~ 10 3 are the Lorentz-factor and the magnetization 
parameter of the wind near the light cylinder. The wind is 
assumed to be monoenergetic at a given latitude. It follows 
from equation (|9h that 



< 7« 



>= 3 7 " 



(10) 



In the particular ca se of the axisymmetric rotator (Bogov- 
alov fc|rCotov 1992|), 



'Jmax 



eB Q R, 
2mc 2 \ 



{RM/c) 



(11) 



2.3 The geometry of the wind flow 



The geometry of the flow and of the IC process are drawn 
in fig. [[]. The neutron star is placed on the axis of rotation 
and ejects the wind radially. The dash-dotted vertical lines 
show the light cylinder. The wind is accelerated in the ac- 
celeration zone by some unspecified mechanism and at the 
distance Rw it has characteristics discussed above. The ac- 
celeration is completed at Rw- Beyond Rw particles in the 
wind move along straight lines without further acceleration. 
The equatorial plane of the pulsar is inclined in relation to 



the observer at the angle a = 33° ( [Hester et al. 1995| ). IC 



phot on s move along the direction of motion of electrons. 



Therefc re only the particles of the wind directed towards 
Earth can produce observable emission. The lines of flow of 
the cold kinetic energy dominated wind after acceleration 
are not exactly radial. Classical mechanics provides a sim- 
ple relationship between the rates of the rotational energy 
losses of the pulsar, E ro t = Ittti, and the angular momen- 
tum losses, L = Jf2: 



E rQ t — L\l. 



(12) 



(13) 



Here I is the momentum of inertia of the neutron star. As 
the electromagnetic field carries off practically no energy in 
the kinetic-energy-dominated wind, it does not carry the an- 
gular momentum either. In the theory of plasma flow in the 
magnetosphere of an axisymmetric ro tator this statement is 
verified immediately ( Bogqyaloy 1997 ) . Under this condition 
it follows from equations (p]),(^), and (p"2|) that after the ac- 
celeration the azimuthal velocity of the wind v v is connected 
with Rw as 

v v _ Rl 
c Rw 

From this relationship and from Fig. |l] it follows that the 
projection of the vector of the velocity of the plasma after 
acceleration lies on the line tangential to the light cylinder. 
The angle 6 between the direction of the motion of rela- 
tivistic particles in the wind and the soft thermal photons 
emitted from the pulsar depends on distance r to the axis 
of rotation: 

sin 6 = cos a— -. (14) 
r 

It is seen from this relationship that at any distance from 
the pulsar there is a non-zero angle between the radial di- 
rection and the velocity of the wind. The Inverse Compton 
scattering of the wind electrons on soft photons emitted as 
a fan-like beam from the inner magnetosphere of the pulsar 
results in the production of hard 7- ray photons. Equations 



source of 
pulsating 
soft emission 



Neutron star 



cone of the 
pulsating soft 




to observer 

Figure 1. A sketch of the trajectories of plasma after accelera- 
tion, and the assumed position of sources of soft photons. 



Q and 



1l3) remain valid for condition j w ^> 7^0 (Bogo- 
valov 1997p 7also fulfilled in the wind with a w ~ 1. Below 
we will not distinguish between the cases with a„, C 1 and 
a w ~ 1 . 

The angular distribution of soft emission close to the 
pulsar is not well known. In the outer gap model the optical 
and soft x-ray emi ssion is generated inside, but n ot far from 
the light cylinder ( Romani fc Yadigaroglu 1995 ). There are 
several sources of soft emission in the magnetosphere corre- 
sponding to the position of the outer gaps. In this model, 
the projection of the motion of photons from these sources 
on the equatorial plane is predominantly directed radially. 
Then the angle between direction of motion of the wind and 
soft non-thermal photons can also be estimated by equation 
([l4|). Note that at Rw 3> Rl the IC flux does not depend 
on the position of the source inside the magnetosphere if the 
soft photons are emitted radially. 

Our calculations of IC radiation are based on the as- 
sumption that the wind is illuminated only by the emission 
that is directed to Earth. However, as we can not exclude 
the existence of an additional component of optical emission 
not directed to the observer, but illuminating the wind, our 
estimates of 7 -ray flux could be considered only as a lower 
limit. 



3 INVERSE COMPTON RADIATION OF THE 
WIND 

3.1 The fields of soft photons 

The soft emission from the Crab pulsar consists of two, 
pulsed and unpulsed, components. The recent ROSAT ob- 
servations allowed one to distinguish the contribution of the 
pulsar in the observed unpulsed flux from the emission of 
the nebula. The main contribution to the unpulsed radia- 
tion of the pulsar comes, most probably, from the thermal 
emission of the hot surface of the neutron star. According 
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Figure 2. The averaged spectral density of photons near the light 
cylinder of the Crab pulsar 



i i nc in!9 -1.11 
+ 1.25 • 10 e exp 



5 • 10" 3 keV 



ph 



cm 3 keV 



(17) 



The photon density of thermal and non-thermal com- 
ponents of low-frequency radiation on the light cylinder is 
shown in Fig. ^. We assume that the density of the both 
photon fields fall off with distance to the pulsar center 7? as 
1/R 2 . 



3.2 The fluxes and spectra of gamma radiation 

The optical depth characterizing the Compton scattering of 
the wind electrons propagating through the radiation fields 
of the pulsar is defined as 



a(uo, 7, v, c)(l — vc)n p h(tj, c, r)dldcu, 



(18) 



where a(ui,~f, v, c) is the total invariant cross-section of the 
IC scattering, uj is the energy of soft photons in mc 2 . In- 
tegration of equation ( |l8| ) over dl is performed along the 
trajectory of a wind electron (see Fig. 1) from the wind for- 
mation position to infinity. 

The 7-ray energy flux at the angle a to the plane of the 
pulsar equator is equal to 



to the ROSAT observations, the total energy flux of the 
unpulsed emission in the range 0.1 - 2.4 keV is 10 34 erg/s 
(Becker & Triimper 1997). The latter could be approximated 
by black-body spectrum with a temperature 1.9 • 10 6 if and 
total luminosity 10 34 erg/s. 

The pulsed soft radiation of the pulsar is dominated by 
nonthermal processes in the magnetosphere. The photons 
in the optical to X-ray band of this radiation play the most 
important role in the production of inverse Compton 7-rays. 
Measurements of th e spectra of the pul s ed optical emission 



1996 



Oke 1969| ; |Perciva] 



by diff eren t groups (Nasuti et al 
et.al. 1993) give the following spectra of optical photons in 
the range 1680 - 7400 A 



F v =3.1( — )" 



mjy, 



(15) 



where v = 6.82 ■ 10 14 Hz. 

The X-ray data in the range of 0.1 - 2.4 keV obtained by 
ROSAT can be approximated by a power-law with the pho- 
ton index 1.5 and luminosity (assuming that the radiation 



is e mitt ed isotropically) 7.1 • 10 erg/s (Becker & Triimper 
19970 

The extrapolation of the soft x-ray emission sp ectrum 
to t he op tical range gives lower flux than is observed (Knight 
1984). This means that the optical emission has a cutoff in 
the ultraviolet region. Unfortunately there are no data for 
the pulsating emission in this region. Therefore we assume 
an exponential cutoff in the spectrum above 0.05 keV. The 
emission is also strongly s uppressed at infrared wavelengths 
(Middleditch et al. 1983). Below we use the following ap- 



proximation of the averaged density of non-thermal photons 
near the light cylinder 



n(e) = Z(t) exp 



8- 10" 



keV Y 



photon 
cm 3 keV ' 



(16) 



where e is the energy of photons in keV and 
Z(e) = 1.82 • 10 16 e~ l f5 + 



N 



4ttD 2 



da 



£ 7 — — (7w(q), uj, E-y, e)n ph (uj, c, r)dwdl, (19) 



where D 

^-( 7w (C),^,£ 7 



2 kpc is the distance to the source, 
is the differential cross-section of the 
Compton scattering of a photon with energy uj and electron 
with 7u, encountering under angle 8 and producing a photon 
with energy E~,. This cross-section wa s obtained after inte- 
gration of the diffe rential cross-section (Berestetskii, Lifshitz 
& Pitaevskii 1971) over the emission angle 

_„2 -1 



da . 
dK^'^ 



21 



(pk) 2 (pk) 2 



+ 



I 3 (I- US 13) 

(■JcUJc) 2 



+ 2{ W) 



{pk)j c u c 
I 



) + 



/ (pk)I 7 c u e (l- US /3) m c 



(pk) 



TVc 



(20) 



where 



(pk) — ~fuj(l — v cos 
n 7 + ^ 



V B = 



7c 



v/l + 2(pfc)' 
1 + (pk) 



-\J a; 2 + C7U) 2 +2~fVL,J cos((?) 



u 



x /l+2(pfe)' 

(Pk) 
l + (pk) ' 



uv \ r 7c 



V / H2p) 

Vc 1 U>c 

Function / in equation ( |2o| ) is defined as follows 

I = (1 - U 2 - 2*7(1 - U)f35 + U 2 (S - f3) 2 y 1/2 . (21) 

Integration of equation ( p"9| ) is performed along the 
trajectory of a electron and over the spectra of soft pho- 
tons. It is assumed that the observer detects photons from 
a monoenergetic beam of electrons from the wind moving 
to the observer. We take 7 max = %E I0 t/N, with 7w(a) = 
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Figure 3. The optical depth of the thermal photon field for the 
IC scattering of electrons with Lorentz factor 7 W = 4 TO 6 (dashed 
line) and 4 ■ 10 7 (solid line). 



7max cos 2 (33°), since the plane of the p ulsar equator is in - 
clined to the observer at the angle 33° (Hester et al. 1995) 



3.3 IC photons from thermal isotropic radiation 

The optical depth r , and therefore the spectra of IC pho- 
tons, depends strongly on the distance R w at which the ki- 
netic energy dominated wind is formed. 

The dependence of the optical depth on i? w is shown in 
Fig. |^ for two different values of the Lorentz factors of the 
wind 7w = 4 • 10 6 and 4 ■ 10 7 Strong dependence of r on 
7w is explained by the fact that the optical depth is dom- 
inated by IC scattering at small distances from the pulsar 
where 2ury(l — cos (#)) S> 1, i.e. the Compton scattering 
takes place in the Klein-Nishina regime. In this regime the 
IC cross-section decreases with the electron energy as j^ 1 , 
which lead to larger optical depth for smaller values of the 
Lorentz-factor of the wind j w . This effect is especially strong 
for small values of R w . In this case the optical depth is ac- 
cumulated by IC scattering close to the light cylinder where 
the collision angle 9 is large, and the IC scattering takes 
place in deep Klein-Nishina regime. 

In many "standard" astrophysical situations the effi- 
ciency of the 7-ray production in the Klein-Nishina regime 
is significantly suppressed because of the synchrotron cool- 
ing of electrons. In the case of a cold wind we deal with 
a unique situation when the synchrotron losses of elec- 
trons are completely suppressed, and thus the wind electrons 
lose their energy only through the inverse Compton radia- 
tion. Note that in deep Klein-Nishina regime, namely when 
2 r y w ui{l — cos (9) > 10 4 , the triplet pair pro duction dom- 
inate s over the inverse Compton scattering ( Mastichiadis 



1991 



Dermer fe Schlikeiser 1991). This process results in 



Figure 4. The spectra of IC radiation of the wind illuminated 
by the thermal emission of the pulsar. It is assumed that the 
kinetic-energy-dominated wind is formed at distances of 1 (solid), 
5 (dashed), and 10 (dotted) light cylinders from the pulsar. The 
curves 1,2,3 correspond to -y max = 1.2 ■ 10 5 ,3 ■ 10 6 ,3 ■ 10 7 . The 
range of observed fluxes of 7-rays in the region above 500 GeV de- 
tected by CAT, CANGAROO, HEGRA and Whipple telescopes 
is shown by a shadowed region. The points with error bars below 
10 GeV correspond to the unpulsed fluxes measured by EGRET. 



we ignore this process here. The same is true also for an- 
other process connected with absorption of TeV 7 -rays in 
the magnetic field of the wind (see Appendix). 

In fig. [I we present the expected 7-ray fluxes of the 
wind. Solid, dashed and dotted lines correspond to the 
fluxes produced by a wind originated at Rw/Rl=1, 5 and 10. 
Remarkably, the spectrum of IC 7 -rays has a specific line 
feature since the radiation is produced in the Klein-Nishina 
regime. Because of this feature, the IC radiation of the wind 
can be easily distinguished from the smooth spectra of the 
Crab Nebula. The corresponding integral fluxes are shown 
in Fig. ^ for different j w and R w . 

The comparison of the calculated spectra with the ob- 
served TeV 7-ray fluxes of the Crab Nebula (for review see 
Weekes et al. 1997) leads to an interesting conclusion that 
for any reasonable Lorentz-factor of the wind , the calculated 
7-ray fluxes significantly exceed the observed fluxes unless 
the wind is formed well beyond the light cylinder. This im- 
plies a meaningful constraint on the 'birthplace' of the ki- 
netic energy dominated wind with ~/ m > 10 6 : R w > 57?l- For 
7 < 1.2- 10 5 , a sim ilar conclusion is imposed by the EGRET 
(Nolan et al. 1993) observations of unpulsed radiation above 
1 GEV. The lack of measurements in the energy region be- 
tween 10 GeV and 300 GeV does not completely exclude a 
possibility of formation of the particle dominated wind close 



production of new electrons. However, as in all interesting 
cases the optical depth r < 1, the secondary electrons do 
not increase significantly the density of the wind. Therefore 



to the light cylinder, if one assumes that 10 < jmax < 10 . 
Although the analysis of the observed synchrotron and IC 
components of the nonthermal radiation of the Crab Neb- 
ula gives certain preference to larger wind Lorentz factors, 
7w ~ 10 6 , it is important to have observational constraints 
on 7? w also for 7 W in the region between 10 5 to 10 6 , which is 
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Figure 5. The integral fluxes of IC radiation of the wind illumi- 
nated by the thermal emission of the pulsar. It is assumed that 
the kinetic-energy-dominated wind is formed at distances of 1 
(solid line ) and 5 (dashed line ) light cylinders from the pulsar. 
The curves 1 and 2 correspond to ^max = 3 • 10 6 and 3 • 10 7 
respectively. The range of observed fluxes of 7-rays in the region 
above 500 GeV is shown by a shadowed region. 



presently missing This should be provided in the near future 
by measurements of the Crab spectrum between 10 and 300 
GeV by the new generation of low-threshold atmospheric 
Cherenkov detectors as well as by the GLAST. 



The source of 
amma-rays 




The spiral formed by 
the narrow beam of soft photons 



Figure 6. The angle between the wavevector of soft photon and 
the velocity of particles in the wind depends on the position of 
the source of soft photons inside the magnetosphere and the di- 
rection of the beam. The minimum efficiency of the gamma-ray 
production is achieved when the source is corotating with the 
pulsar. The narrow beam of non-thermal soft photons produces 
a spiral in space. The point of crossing of this spiral with the 
trajectory of an electron moving to an observer is the source of 
7-ray photons, which moves with a velocity very close to the light 
velocity. At small angles between the beam of the soft photons 
and the velocity of particles in the wind, the delay in the arrival 
time of hard photons and optical photons is small. The width 
of the gamma-ray pulse is also small. The phase curves for opti- 
cal radiation illuminating the wind and in VHE gamma-rays are 
expected to be similar. 



3.4 The interaction of the wind with the 
nonthermal emission 

The IC 7 radiation caused by illumination of the wind by 
pulsed soft photon emission should be modulated at the 
same pulsation period. This conclusion is true even for the 
isotropic wind. For calculation of the IC fluxes it is assumed 
that the nonthermal source of soft emission is located in- 
side the light cylinder (see Fig. Q). It is easy to show that 
under realistic conditions the phase curve of the gamma- 
rays should be close to the phase curve of the soft emission. 
Indeed, let us assume for simplicity that the beam is very 
narrow (delta- function like). Owing to the rotation of the 
pulsar, it will form a spiral. The point where this spiral 
crosses the line tangential to the light cylinder and directed 
towards the observer is the source of 7-rays (see Fig. ^|) . The 
width of the 7-ray pulse is determined by the delay in the 
arrival time St of optical and gamma-ray emission to the 
observer: 



where 9 is the maximal value of the angle between the wave 
vector of photons and the velocity of the particles in the 
wind, T is the period of rotation. In the most interesting 
case, 8< 1, the phase curves of gamma-ray emission pro- 
duced in the process under consideration and those of the 
optical emission illuminating the wind are expected to be 



similar. However, because we observe only the pulsed soft 
emission directed to us and do not observe other possible 
part of the emission (not directed towards the observer) that 
can illuminate the wind, there could be a difference between 
the 7-ray and directly observed soft emission light curves. 

It is seen from Fig. ^| that the average density of non- 
thermal photons is several orders of magnitude larger than 
the density of the thermal photons. Therefore, the IC opti- 
cal depth for an electron is much larger than in the case of 
thermal radiation (see FigJ^). In particular, at Rw < 5Rl 
the optical depth r > 1. As the IC scattering takes place in 
the Klein-Nishina regime, one or two interactions are suffi- 
cient to destroy the wind, and the whole energy of the wind 
electrons would be transferred to 7 ray emission with huge 
luminosity > 10 38 erg/s. This very fact excludes the pos- 
sibility of formation of the kinetic dominated wind within 
5Rl- Moreover, this conclusion can be extended to larger 
distances. Fig. ^ shows the spectra of the emission for two 
maximal energies of the electrons in the wind. Solid lines 
show the spectra for ^ m ax = 3 • 10 6 and dashed lines show 
the spectra for "f ma x = 3 ■ 10 7 . In contrast to the IC 7- 
rays from interaction with thermal emission, there are no 
lines in the spectra of the emission because the IC scatter- 
ing takes place in the Thomson regime. The spectra consist 
of two well separated broad components, corresponding to 
two components in the spectra of the soft nonthermal emis- 
sion presented in Fig. 0. The comparison of calculated 7-ray 
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Figure 7. The optical depth for IC scattering of electrons on 
nonthermal pulsating radiation of the pulsar. It is assumed that 
the nonthermal photons are directed radially. Curves are given 
for 7u, = 4 • 10 6 (dashed line) and for -f w = 4 ■ 10 7 (solid line). 



fluxes wit h upper limits on pu lsed emission re ported by the 
Whipple ( |Weekes et al. 1998[ ) and HEGRA ( Aharonian et 
al. 199i|) groups shows that the wind must be formed at 
relatively large distances from the pulsar, Rw > 30Rl- 

The spectra shown in Fig. |^ were calculated at the as- 
sumption that the nonthermal source emits photons along 
the radial dir ection from the pulsar as it hap pens in the outer 
gap model (Romani & Yadigaroglu 1995). Self-consistent 



MHD solutions confirm that the plasma should move pre- 
dom inan tly radially in the pulsar magnetosphere ( Bogovalov 
199?|)Hrherefore the source of nonthermal photons should 
indeed emit along the radial direction from the pulsar. In the 
model by Lyubarskii (1996) the source of the soft nonther- 
mal photons emits in the direction opposite to the direction 
of rotation of the pulsar. If so, the 7 -ray emission would be 
even higher than shown in Fig. ^. 

We are not aware of models in which the source of pho- 
tons corotates with the pulsar and emits photons tangen- 
tially along the direction of rotation. In this case the angle 
of interaction 8 is minimum and correspondently the IC flux 
is expected to be reduced. Although it is almost unlikely that 
the plasma (and source of soft photons ) can corotate with 
the pulsar, it is perhaps wise to consider this limiting case 
from pedagogical point of view in some detail. 

Optical and soft X-ray emission can be produced by 
plasma moving with a Lorentz-factor of order 7 e ~ 10 2 — 10 3 
in the pulsar magnetosphere. At corotation with the pulsar 
the plasma have azimuthal velocity v v = — , where r is the 
radius in the cylindrical system of coordinates. Since these 
particles should also be directed towards the observer, they 
have the component of the velocity v z = ccosa. From the 
relativistic relationship it follows that 



Figure 8. The spectra of pulsating TeV gamma-ray emission 
produced by the wind with ^/ m ax = 3 ■ 10 6 (solid lines) and 
7 moI = 3 ■ 10 7 (dashed lines) Numbers 1,2 and 3 near the solid 
lines correspond to R w /Rl = 30, 50 and 70, respectively. The 
curves 4 and 5 correspond to R w /Rx J = 70 and 100. The upper 
limits on the pulsed radia tion from Crab Neb ula above 300 GeV 
given by Whipple group (Weeke s et al. 1998) and abo ve 1 TeV 



given by HEGRA collaboration (Aharonian et al. 1999) are also 



shown. The beam of soft photons is assumed to be radial. 



where r B is the distance from the source to the center of 
the pulsar, v r is the radial component of the velocity. The 
soft photons are emitted at the angle ip as shown in Fig. [| 
Neglecting the term I/7 2 we obtain that 



sin 7/) 



(24) 



7, 



= 1 + 7e[cos a + ( ) +( — )], 



(23) 



This relationship shows that a source that emits soft photons 
in direction to the observer and corotating with the pulsar 
can not be located exactly on the light cylinder. 

Fig. |^ demonstrates how the spectra of IC radiation 
depend on the position of the corotating source of the soft 
photons inside the pulsar magnetosphere for parameters of 
the kinetic energy dominated wind at R w = AQRl- Curve 
1 corresponds to the emission from the corotating source 
located on the surface of the pulsar. The same flux is pro- 
duced by the source of soft photons with the beam directed 
radially from the pulsar; the IC flux generated from the last 
source does not depend on r s . Therefore comparison of the 
curve 1 with others allows one to compare the fluxes of IC 
photons produced at the scattering of electrons on the soft 
emission from the corotating source and the source with the 
soft photons emitted radially, but located on the same dis- 
tance r s . It follows from this figure that the flux of the IC 
photons can be reduced to zero only at the position of the 
source of soft photons not far from the light cylinder. At 
any other position of the source the efficiency of generation 
remains high at r s < O.GRl- The existing upper limit on the 
pulsating flux from the Crab pulsar means that either the 
wind is formed beyond 30-Rz, or the wind is formed close to 
the light cylinder but the pulsating source of soft nonther- 
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Figure 9. The dependence of the differential spectra of 7 -ray 
emission on the position of the source of soft photons inside the 
pulsar magnetosphere for parameters of the wind •y-max = 3 • 10 6 
at Rw = 40Rl- The source of soft photons is placed on the 
distance r a = 0.,0.3, 0.5, 0.6Rl respectively for curves 1, 2, 3, 4. 
The direction of the beam of the soft photons was defined under 
the assumption of corotation of emitting particles with the pulsar. 



mal photons corotates with the pulsar and is located close 
to the light cylinder. 



4 SUMMARY 

The Crab Nebula is a unique cosmic laboratory with an 
unprecedentedly broad spectrum of the observed nonther- 
mal radiation that extends throughout 21 decades of fre- 
quency - from radio wavelengths to very high energy 7-rays 
(for review see e.g. Aharonian & Atoyan 1998) It is com- 
monly accepted that the synchrotron nebula is powered by 
the relativistic wind of electrons generated at the pulsar and 
terminated by a standing reve rse shock wave at a distance 
r s ~ 0.1 pc ( [Rees fc Gunn 1974| ) . The relativistic MHD mod- 
els, evep in their simplified form (e.g. ignoring the axisym 



associated with the nonthermal and thermal low-energy ra- 
diation components of the pulsar, respectively. 

The unpulsed component of 7-ray emission associated 
with thermal radiation of the pulsar with temperature ~ 
2 • 10 6 K is produced in deep the Klein-Nishina regime, and 
therefore has a very sharp (line-like) spectral feature which 
peaks at energy E ~ j w m c c 2 . Detection of this compo- 
nent would therefore result in unique information about the 
Lorentz-factor of the bulk motion of the wind. The nonther- 
mal radiation of the pulsar has rather broad energy spec- 
trum which extends to optical and infrared wavelengths, and 
therefore the IC 7-ray emission associated with this compo- 
nent takes place, to a large extent, in the Thomson regime. 
This results in a broad 7-ray spectrum with a sharp cutoff 
at E ~ 7™ m c c 2 . 

The absolute 7-ray fluxes of both components depend 
strongly on the site of formation of the kinetic dominated 
wind, as well as the Lorentz-factor and the geometry of prop- 
agation of the wind. Thus even the flux upper limits of these 
7-ray components should provide important constraints on 
the wind parameters. In particular, we show that the com- 
parison of the calculated flux of the unpulsed inverse IC 
emission with the measured 7-ray flux of the Crab Nebula 
excludes the possibility of formation of the kinetic-energy- 
dominated wind within 5 light cylinder radii of the pulsar, 
i? w > 5Rl- The analysis of the pulsed IC emission, calcu- 
lated under reasonable assumptions concerning the produc- 
tion site and angular distribution of the optical pulsed radi- 
ation, yields even tighter restrictions, namely J? w > 30-Rl- 

The mechanism of 7-radiation of the wind of the Crab 
pulsar discussed in this paper should certainly take place 
in other pulsars as well. However, from the point of view 
of detection of this radiation, the Crab is a unique object 
due to its very powerful wind and relatively high luminosity 
of thermal and nonthermal low-energy radiation, which pro- 
vides seed photons for the IC scattering. In other pulsars the 
IC 7-ray fluxes of unshocked winds are expected to be below 
the detection threshold of current 7-ray instruments, unless 
the kinetic energy dominated winds of pulsars are produced 
very close to the light cylinder. The situation could be dif- 
ferent in binary systems containing a pulsar and luminous 
optical companion, the latter being an effective supplier of 
seed photons for IC scattering. For example, the pulsar/Be 
star binary system PSR 1259-63 seems to be a unique ob- 
ject for the search for IC TeV radiation from b oth shocked ( 



metric structure of the wind and its interaction with "the Kirk Ball fc Skjaeraasen 199S| ) and unshocked ( J Kirk fc Ball 



optical Ifilaments 1 ) . successfully describe the general charac- 
teristics of the synchrotron nebula, and predict realistic dis- 
tributions of relativistic electrons and magnetic field in the 



down stream region behind the shock (Kennel & Coroniti 



1984). Meanwhile our knowledge of the unshocked wind, i.e. 
about the region between the pulsar magnetosphere and the 
shock is based only on theoretical speculations. Moreover, 
it is commonly thought that the wind could not be visible 
in the region upstream of the termination shock because 
the relativistic electrons and magnetic field in wind move 
together, thus the unshocked wind does not produce syn- 
chrotron radiation. In this paper we show that the kinetic 
energy dominated wind nevertheless could be directly ob- 
served through its IC radiation because of the illumination 
of the wind by low-energy radiation of the pulsar. The 7-ray 
emission consists of two components, pulsed and unpulsed, 



199E) winds of the pulsar. 
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APPENDIX A: OPACITY OF THE MAGNETIC 
FIELD OF THE WIND FOR THE TEV 
PHOTONS. 



In the paper of (Cheung & Cheng 1994) it was argued that 
TeV photons produced near the light cylinder of the Crab 
pulsar will be absorbed, because of conversion of these pho- 
tons in e pairs in nonuniform magnetic field. We show here 
that accurate estimate of the rate of the conversion of pho- 
tons in pairs taking into account electric field existing in 
the wind give negligible absorption of TeV gamma-rays pro- 
duced via IC scattering. 

The motion of the wind occurs under frozen in condi- 
tion. Particles move not only along field lines, but there is 



Direction 
of emission 
of ICS photon 




Figure Al. The geometry of the MHD plasma flow. 



e A 




B 

Figure A2. The geometry used in equation ( |A2| ). 



component of the velocity Vd directed perpendicular to the 
magnetic field line. This is so called drift velocity defined by 
the expression 



E x B 

— w 



(Al) 



Vd is comparable with c beyond the light cylinder, where 
the electric and magnetic fields are of the same order. The 
schematic relationship between the magnetic and electric 
fields and the velocity of plasma is shown in Fig. Al. In this 



situation an energetic photon is emitted along the velocity 
vector of the particle at large angle to the magnetic field. 
However, this happens in the region of crossed electric and 
magnetic fields, where the coefficient of absorption of the 
photon is modified by the electric field. The probability of 
conversion of the photon on unit length, taking into ac coun t 
the electric field, was defined by Daugherty & Lerche: (197 



£ = 0.23a 



mc B X (1-E 2 /B 2 ) 
~~h~B^ (1- Er^/B) 



exp 



8 B c 



3 E^xB 



where 



B' 



EL 

B 2 



(A2) 



(A3) 



and a — 1/137. r/ x and % are the components of the unit 
vector directed along the velocity of the photon as it is shown 
in Fig. and E 1 is the energy of photon expressed in units 
of mc 2 . 

The components of the velocity of the particle in the 
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system of coordinates presented in Fig. A2 are as follows the wind are not converted into pairs even at energies ~ 100 



Vx = Vy = 0, v z = c^Jl - §£. The corresponding com- TeV and CTw ~ L 
ponents of the unit vector r\ of the photon at the place of 
the emission will be 77* = 2J=. = £2, r\ y = 0, ^ = y/l - rj£. 
Substitution of these components in equation ( |A3| ) gives us 

_ E c 

X ~ B-y 2 (l + v/c)v [ ' 

This factor is extremely small, \ ~ 10~ 14 , for the expected 
parameters of the wind . 

Photons propagate along straight lines. Charged par- 
ticles move in electromagnetic fields and their trajectories 
diverge from straight lines. If photon would moved along 
the trajectory of the charged particles, they would never 
convert in pairs. Therefore it is clear that the probability 
of conversion of a photon into a pair basically depends on 
the radius of curvature R c of the trajectory of the charged 
particles, but not the magnetic field. 

The angle $ between the direction of propagation of the 
photon and direction of propagation of the emitting electron 
depends on the path length of the photon I; 

*=■!. (A5) 

The dependence of components of the vector r\ on 9 is 
as follows 

r/x = r/xo cos 1? + r/ z o sin # sin ip (A6) 

rj v = sin ■& cos tp (A7) 

Neglecting by terms of the order I/7 2 and sin" 1? in powers 
higher than 2, we obtain the following estimate for x at small 



*~sV'-(f)- <A8) 

As we are intereste d in the maximal values of x> we can 
neglect the term yl — (-§) 2 in this expression. Assuming 
that after the light cylinder the total magnetic field decreases 
as 1/r, we obtain for the function q — 3 g B g x the following 
estimate of the upper limit 

8B CI R C . . 

1 = T^r-g , (A9) 

and for the probability of conversion of the photon into a 
pair we obtain finally the estimate 

t/ -„„ rnc Bi c r L 8B CI R C 

* * °' 23Q - b-rZ exp ( " 3iM^ (A10) 

where B\ c is the magnetic field on the light cylinder and R c 
is also taken on the light cylinder. 

The solution of the problem of the relativistic wind flows 
in the model of an axisymmetrically rotating star allows us 
to estimate the radius of curvature R c of the trajectories of 



the particles in the wind (Bogovalov 1997) as 

■i- = ^r- ■ (All) 

R c 7 2 r v ; 

This estimate shows that the radius of curvature of the 
trajectories of particles is much more than the curvature ra- 
dius of the magnetic field line, which is of the order r. Such 
a large curvature radius gives 1 for almost any param- 
eters of the wind. This is why the IC photons produced in 
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